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ABSTRACT 

Leading QCD vacuum polarization contributions to the electroweak parame- 
ter Ap are evaluated numerically using several different prescriptions for the 
gluon self-energy. Simple theoretical estimates of the asymptotic behavior are 
given. The results show a significant contribution from the leading infrared 
renormalons when Ap is expressed in terms of the top-quark pole mass and its 
absence when the MS running mass is employed. The calculations are applied 
to estimate higher order QCD contributions to Ap. 
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The higher order effects arising in the ratio M/rh(M) (M is the pole mass of a quark and 
rh the MS running mass) from a single chain of vacuum polarization diagrams has been 
recently studied in considerable detail |], |j. In particular, in Ref. |] it has been empha- 
sized that the basic building block in these calculations is a gauge-dependent amplitude, 
and therefore there is a considerable latitude in its definition. This can be explicitly illus- 
trated by considering the one-loop gluon vacuum polarization diagram evaluated in the £ Q 
Background-Field-Gauge (BFG) 0: 
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where b = 11 — 2n//3 is the coefficient of the leading term in the QCD /3-function, n/ the 
number of massless fermions, k the external gluon momentum, £ Q is the gauge parameter 
associated with the quantum loops, and the superscript MS means that the MS renormal- 
ization has been carried out. An important feature of Eq.(|l]) is that it includes correctly 
the logarithms associated with the running of a s . However, the accompanying constant 
depends on the choice of £ Q . Most current evaluations of QCD vacuum polarization effects 
are based on the leading rif approximation, coupled with a non-abelianization prescription 
0: starting from the complete fermionic contribution to the gluon self-energy, one includes 
the dominant gluonic component by replacing —2nf/3 — > b. This is equivalent to the ne- 
glect of the 1/6 term in the expression between curly brackets in Eq.(^J). For brevity, we 
will refer to this procedure as the "non-abelianization" (NA) approach. In Ref. [|l| several 
different prescriptions for the basic building block were discussed. In this paper we extend 
these considerations to evaluate the higher-order effects induced by a chain of vacuum po- 
larization graphs in the QCD corrections to (Ap)f, the fermionic component of Ap. We 
recall that this amplitude plays a significant role in electroweak physics, as it contains the 
leading contributions, for large M t , of the basic corrections Ar, Af, and Ap H, [IJ. 

A very convenient framework to carry out this analysis is provided by the work of Smith 
and Voloshin (S-V) 0. Setting M b = and writing 

Of 

(A P ) f = ^- 2 M?(l + 5 QCD ), (2) 

where M t is the top-quark pole mass and 5 QCD is the QCD correction, these authors have 
shown that the one-loop vacuum polarization insertion in the 0(a s ) correction is given by 

(<WW = f^ r[w(x)+s(x)](l + U(-K 2 ))dx 2 , (3) 
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where x = k/p, k 2 = —k 2 is the invariant euclidean gluon momentum, and the weight 
functions w(x) and s(x) are given by 
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The function w(x) occurs in the combination of vacuum polarization contributions to 
(Ap)/ involving the (t,b) isodoublet, while s(x) is associated with the mass counterterm 
corresponding to the M t pole definition. Defining the leading vacuum polarization effects 
as the contributions of 0(a™) involving n — 1 one-loop vacuum polarization insertions 
(i.e. the maximum number of insertions), it is clear that an analogous approach can be 
employed in their evaluation. As the integral in Eq.(^) is convergent, it suffices to replace 
1 + k 2 ) — > (1 — n(— K 2 )) -1 with rif = 5 (i.e. b = 23/3), expand the geometric series 
and carry out the integrations numerically. 
Using the generic expression: 



where a = a s /n, the coefficients p n are displayed in Table 1 up to n = 10 for five different 
prescriptions to evaluate the basic building block U(k 2 ): 1) NA, in which case the p n coef- 
ficients are denoted as a n ; ii) Pinch Technique (PT) [[/]] which leads to the same expression 
as the £ Q = 1 BFG (p n — > p n ); iii) £ Q = (or Landau BFG) with p n — ► £ n ; iv) £ Q = —3 
BFG, which corresponds to the minimum of the non-logarithmic contribution in Eq.(|T|) 
(p n — > v n ); v) MS running coupling a s (n) (p n — > r n ). It was pointed out in Ref.|l| that, 
by a rather remarkable coincidence, the £ Q = —3 self-energy equals the 0(a s ) correction 
to the QCD potential between two infinitely heavy quarks ||, and it is thus amenable 
to a physical interpretation. We have also included the elementary running-coupling pre- 
scription, in which only the logarithmic term in Eq.(^) is retained. Although it cannot be 
obtained from Eq.(|I]) for any real value of £ Q , this approach is perhaps the most intuitively 
simple, and it is frequently employed in the literature || [10| . Another attractive feature 
is that a s (n) is gauge- invariant. 

The evaluation procedure involves a numerical integration up to a cutoff A, with the 
tail beyond the cutoff evaluated analytically using the asymptotic expansion for the weight 
functions w(x) and s(x). Different values of A were employed to check the stability of 
the evaluation, and a similar method was employed for the very small x integration re- 
gion. Typically, the variation found is of the order of a few parts in 10 8 . An important 
consistency check involves the a\ coefficient, as this must lead exactly to the 0(a 2 (M t )) 
contribution proportional to n/. Replacing b — > — 2n^/3 in the n = 1 term of Eq.(^j), 
we have (8/9)ain/a 2 (M 4 ). Keeping eight significant figures our numerical evaluation gives 
a\ = 2.0094894, which leads to 1.7862128 nfa 2 (M t ). To this accuracy, this value coincides 
with the very precise determination reported in Ref.[|ll|]. Other consistency checks involve 
the predicted values of the asymptotic coefficients, which will be discussed later. 

A glance at Table 1 shows that all the coefficients are positive and the various sequences 
converge to well-defined limits, with the asymptotic behavior setting in remarkably quickly. 
This is analogous to the result for M t /m t (M t ), obtained by analytic methods P], 0|. The 





structure of Eq.(|6]), involving the factorial enhancement, the asymptotically constant p n 
of fixed sign, and the expansion variable ba(M t )/2, shows the characteristic Borel non- 
summable behavior associated with the leading infrared-renormalon contribution. In our 
case, this is the pole in the Borel plane at u = btj (4-7r) = 1/2, where t is the Borel parameter 
[T], 0. Physically, this behavior emerges because in the evaluation of Feynman amplitudes 
the gluon self-energy (or alternatively the strong interaction coupling) becomes very large 
at low values of k/p, and perturbation theory breaks down fT0|j . 

Combining the above results with those obtained in Ref.flJ], one can readily find the 
corresponding expansions for A QCD , the QCD correction when Ap is expressed in terms of 
mf(M t ): 

Of 

(Ap) f = ^J^m 2 (M t )(l + A QCD ). (7) 

Defining again the leading vacuum polarization effects as the contributions arising, in a 
given order in a s , from the maximum number of vacuum polarization insertions, we have 

M t V 8 ~ (ba(M t )\ n , 

1 1 + - fl (M t ) £ A*J — n\, (8) 
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where the p n coefficients are evaluated, for the i)-iv) prescriptions, in Ref. [[[[]. We have 
also calculated the coefficients p n for the running coupling approach, by the same analytic 
methods. Combining Eqs.(|l]) and one obtains for the leading contributions 



Aqcd = ~-a(M t ) £ Pn 1 n\, (9) 

with p n = p n - Pn- Using the notation p n — >■ pn,a n J n ,v n ,r n for the PT, NA, £ Q = 0, 

^ Q = —3, and running coupling prescriptions, respectively, these coefficients are displayed 

in Table 2 up to n = 10. Comparing Tables 1 and 2, we see that the pattern of the M t 2 and 

rht(Mt) expansions is dramatically different. In the latter case, the coefficients are much 

smaller, decrease as n becomes large and the signs alternate after n = 2. 

The general features and many of the detailed results displayed in Tables 1 and 2 can 

be understood on the basis of very simple asymptotic estimates. Expressing the gluon 

self-energy in the abbreviated form II(— n 2 ) = (a(M t )b/A) \n(l/( 2 x 2 ), where ( 2 = e~ p and 
p2 = 5/3 + ( 4/ 6)[1 _ 3(1 _ £ Q ) (7 + £ q)/16]) we have 



(<W W = \<M t ) f] J™ 2x [w(x) + s(x)) In" (1) dx. (10) 



n=0 

As the infrared-renormalon contributions involve the very small x region, we introduce an 
ultraviolet cutoff x max and employ the expansion 2x[w(x) + s(x)} = —4 + |x 2 + 2x 3 + . . . 
{x <^ 1). Setting y = \n(l/(x) as integration variable, and choosing x max = 1/C (where 
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n(— n 2 ) = 0), the integrals reduce to factorials and we obtain 
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These contributions are associated with poles on the positive real axis of the Borel plane 
(infrared renormalons) . The first term arises from the pole nearest to the origin (u = 
1/2) and gives the leading asymptotic behavior for large n, while the others involve more 
distant renormalons (u = 3/2, 2 . . .) and represent subleading contributions. To study the 
dominant asymptotic contributions from the ultraviolet region we introduce an infrared 
cutoff x m i n in Eq.flTU) and use the asymptotic expansion 2x[w(x) + s(x)] = 24(lnx — 
l)/x 3 + (157 — 1921nx)/x 5 + . . . (s > 1). Choosing x min = x max = 1/C and proceeding as 
before, we obtain 
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The two terms between curly brackets correspond to the leading and next-to-leading ul- 
traviolet renormalons (poles on the negative real axis of the Borel plane). The unusual 
contributions proportional to n within the curly brackets arise from the additional loga- 
rithms from the w(x) expansion, and correspond to singularities of the form u/(l + u) 2 , 
uj (2 + u) 2 in the Borel plane (i.e. combinations of simple and double poles). Eqs. (|6|, |llT , |T2|) 
imply p n — > 1/C = e p I 2 as n — » oo, the same asymptotic behavior as the fi n 0, so 
that p n — > 0. In fact, it is very instructive to compare Eq . (|TTD with the corresponding 
expansion in M 2 /fn 2 (M). One readily discovers that the two leading infrared-renormalon 



expansions have equal magnitude and opposite sign and therefore cancel in A QCD \ 
bining Eqs. ( fn|jr2] ) with the corresponding expansions for M 2 /m 2 (M), one obtains 
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(A^C^^E(^) n n,^ + ... (14) 

We see that, in contrast with 5 QCD , the leading asymptotic behavior in A QCD is associated 
with the nearest ultraviolet renormalon (u = —1), which explains the alternating signs in 
Table 2 and the much reduced coefficients relative to Table 1. The simple asymptotic for- 
mulae of Eqs. (jnyT3) approximate the p n coefficients (Table 1) with considerable accuracy 
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for n > 2. On the other hand, Eqs. (13,14) give an accurate estimate of the p n coefficients 
(Table 2) for n > N, where iV depends on the prescription employed. We illustrate this 
for the v n and v n coefficients in the last column of Tables 1 and 2. We see that the v n and 
v n fits are quite precise for n > 1 and n > 2, respectively. We have extended this compar- 
ison up to the 17th coefficient and find that the numerical and asymptotic evaluations of 
Vn = 9.79 x 10 -5 differ by only 1.6x 10~ 8 . As the leading asymptotic contribution for A QCD 
arises from the u — — 1 ultraviolet renormalon, it is Borel summable and can be expressed 
in terms of ordinary functions and exponential integrals. Subtracting the n — 0, 1 terms 
from this sum and using the £ Q = —3 prescription, we find that, for M t = 175GeV, it 
contributes only 4.3xl0 -4 to A QCD . Thus, in contrast with the large effects from infrared 
renormalons in 5 QCD , the n > 2 terms in the leading ultraviolet renormalon expansion give 
a small contribution to A QCD . 



In the limit Mb = 0, recent exact calculations [11, 12] lead to the expansion 



S QCD = -2.8599 a{M t ) - 14.594 a 2 (M t ) + ... (15) 

Using Eq.@ and pi, a\, £\, v\, and r\ from Table 1, the 0(a 2 (M t )) coefficients from the 
vacuum polarization contributions are: i) —23.401 (PT), ii) —20.541 (NA), iii) —19.647 
(£ Q = 0), iv) —14.821 (£ Q = —3), v) —11.406 (running coupling). We see that the vacuum 
polarization contributions are of the same sign and roughly the same magnitude as the exact 
calculation. However, the PT, NA, and £ Q = prescriptions overestimate the 0(a 2 (M t )) 
coefficient by 60.3%, 40.7%, and 34.6%, respectively. The fact that the NA does not 
provide an accurate determination is related to the observation that when one applies the 
BLM rescaling Jl3| in the Ap rather large residual 0(a 2 (M t )) term remains [Q. 



On the other hand, the £ Q = —3 prescription approximates this coefficient very accurately, 
within 1.6%. The leading £ Q = —3 calculation also gives an 0(a 2 (M t )) coefficient —8.476 
in rht(M t )/M t and +8.746 in the inverse M t /rht(M t ) expansion. The exact results (in 
the Mb = limit) are —9.125 and +10.903 respectively, so that the £ Q = —3 approach 
underestimates the answer in these cases by 7.1% and 22.3%, respectively. In summary, 
we have seen that the £ Q = —3 BFG presents a number of interesting features: i) within 
the BFG it generates the smallest constant accompanying the logarithmic term, and is 
expected to provide the smallest asymptotic coefficients in infrared-renormalon effects; ii) 
it coincides exactly with the 0(a s ) correction to the potential between two infinitely heavy 
quarks 0; iii) it approximates very well (within 1.6%) the 0(a 2 (M t )) coefficient in 6 QCD ; 
iv) it approximates quite well the corresponding term in the rh t (M t )/M t expansion (within 
7.1%), although it is less precise (22.3%) in the inverse M t /rh t {M t ) case. In contrast with 
the other approaches, the running coupling prescription underestimates the 0(a 2 (M t )) 
coefficient in 5 QCD by ~ 22% (although it fares much better if the comparison is made with 
the contribution to the exact result not involving the double triangle graph, in which case 
it overestimates the answer by 9.1%). 

Assuming that the v n series (Eq.(D) with p n —>■ v n ) provides also the dominant contribu- 
tions for n > 1, we can employ it to estimate the higher order effects in 5 QCD . We note that 



5 



this expansion involves large higher order coefficients: for example the coefficients of a 3 (M t ) 
and a 4 (M 4 ) are 81.3 and 1,435.2, respectively. As an illustration, setting M t = 175GeV 
(c4 5) (175GeV) = 0.10744 for a[ 5 '(m z ) = 0.118) and summing up to the previous to small- 
est (n = 6) or smallest (n = 7) terms, we find — 5 QCD = 0.1219 or 0.1223 (in the evaluation 
we employ the exact 0(a 2 (M t )) coefficient —14.594 rather than —14.821). This is close, 
albeit somewhat higher, than the value 0.1200 found in Ref. |]J by using the BLM opti- 
mization of Mt/rht(Mt), together with other considerations. It is also somewhat larger 
than the values 0.1192, 0.1198, and 0.1197, obtained by applying directly the BLM 
fl3|l , PMS fill ) and FAC [0 optimization procedures to Eq . (|15|) , and larger by 0.0070 



(0.0075) than the value 0.1149 from Eq . fll5|) . As the coupling constant prescription gives 
the smallest renormalon effects in Table 1, we have repeated the above calculation using 
the r n coefficients (retaining again the exact 0(a 2 (M t )) term) and find — S QCD = 0.1195 
(n = 5) and 0.1198 (n = 6), which are even closer to the optimized values mentioned 
above. It is worth noting that the 0.1149 value from Eq.(|l5|) is marginally consistent with 
the 5.3 xlO -3 error estimate around 0.1200 given in Ref.0. The renormalon calculations 
described above favor, however, values of 5 QCD closer to the central value 0.1200 obtained 
in that work. 

In summary, the numerical evaluation and theoretical asymptotic estimates of leading 
vacuum polarization contributions indicate the presence of a significant contribution from 
the dominant infrared renormalons in Ap when it is expressed in terms of the top pole mass, 
and its absence when the MS running mass is employed. The detailed results depend 
on the precise definition of the basic building block, the gluon self-energy. One of the 
vacuum polarization prescriptions, the £ Q = —3 BFG, leads to an 0(a 2 (M t )) coefficient 
that approximates very accurately the result from the exact calculation. Assuming that 
the infrared-renormalon effects evaluated with this prescription also dominate the higher 
order contributions, one finds that they increase the value of S QCD relative to Eq.([15|) by 
an amount that is close to the effects previously obtained by optimization considerations 
[TJ, [Ll| [T7|]. A similar conclusion is reached if one simply employs the running coupling 
olAk) in the evaluation of the renormalon effects. 
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n 


Pn 


a n 






r n 


/ \ 

{V n )as 





1.07247 


1.07247 


1.07247 


1.07247 


1.07247 




1 


2.28926 


2.00949 


1.92206 


1.44994 


1.11577 


1.40672 


2 


2.50023 


1.93952 


1.77927 


1.03706 


0.63733 


1.02925 


3 


2.93381 


2.35629 


2.20476 


1.59249 


1.33432 


1.58413 


4 


2.88349 


2.19666 


2.01084 


1.18702 


0.73426 


1.18785 


5 


3.00202 


2.34266 


2.17124 


1.47723 


1.17998 


1.47672 


6 


2.96229 


2.26905 


2.08517 


1.29537 


0.88504 


1.29543 


7 


2.99567 


2.31705 


2.13967 


1.40676 


1.06946 


1.40673 


8 


2.98016 


2.29113 


2.10957 


1.34168 


0.95921 


1.34169 


9 


2.98995 


2.30623 


2.12695 


1.37910 


1.02336 


1.37910 


10 


2.98490 


2.29800 


2.11740 


1.35804 


0.98684 


1.35804 



Table 1: The coefficients p n in Eq.(6). The p n , a n , £ n , v n , r n correspond to the pinch- 
technique, non-abelianization, £ Q = 0, £ Q = —3, and a s (K) prescriptions for n(— k 2 ), respec- 
tively. The last column gives the result of the asymptotic estimates (Eqs. (12,13)) for the 
£ Q = —3 case. The normalization of Eq.(6) corresponds to p = (1 + 7r 2 /3)/4 = 1.07247... 



n 


Pn 








^'n 


(Vn)as 





0.07247 


0.07247 


0.07247 


0.07247 


0.07247 




1 


-0.49403 


-0.33359 


-0.36905 


-0.20835 


-0.37948 


-0.09436 


2 


-0.34464 


-0.26331 


-0.24266 


-0.19952 


-0.12460 


-0.20072 


3 


0.00058 


0.06757 


0.08718 


0.17430 


0.22649 


0.16495 


4 


-0.06959 


-0.07918 


-0.08550 


-0.14297 


-0.20541 


-0.14241 


5 


0.02159 


0.04066 


0.04737 


0.09628 


0.15011 


0.09578 


6 


-0.01921 


-0.02726 


-0.03084 


-0.06154 


-0.09942 


-0.06148 


7 


0.00924 


0.01521 


0.01757 


0.03712 


0.06176 


0.03710 


8 


-0.00576 


-0.00890 


-0.01023 


-0.02178 


-0.03688 


-0.02178 


9 


0.00308 


0.00496 


0.00574 


0.01247 


0.02141 


0.01247 


10 


-0.00174 


-0.00278 


-0.00321 


-0.00703 


-0.01219 


-0.00703 



Table 2: Same as in Table 1, for the coefficients p n in Eq.(9). The normalization of Eq.(9) 
corresponds to p = p — 1 = 0.07247 . . . 
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